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Abstract 

Using Monte Carlo method we study a two-dimensional model with infinitely many 
absorbing states. Our estimation of the critical exponent /3 ^ 0.273(5) suggests that the 
model belongs to the (1+1) rather than to (2+1) directed-percolation universality class. 
We also show that for a large class of absorbing states the dynamic Monte Carlo method 
leads to spurious dynamical transitions. 

Recently, nonequilibrium phase transitions have been very intensively studied. To some 
extent this is motivated by their various potential applications ranging from catalysis Q to 
epidemic processes Q to interfacial behaviour in random media Q. Another motivation is a 
desire to classify a rich behaviour of these systems into some universality classes in a manner 
resembling a relatively complete classification of equilibrium phase transitions. The basic 
idea of such an approach is to identify a parameter which determines to which universality 
class the model actually belongs to, as e.g., dimensionalities of the order parameter and of 
the embeding space for the equilibrium phase transitions. It was already suggested that for 
nonequilibrium models the corresponding parameter might be the number of absorbing states 
and in particular models with a single absorbing state should belong to the so-called directed- 
percolation (DP) universality class 0]. Up to now there is an ample numerical evidence in 
support of this conjecture [^. 

There are also some indications that certain other models might belong to the so-called 
parity-conserving universality class and most likely models with two absorbing states belong 
to this class |6|, ^. Little is known, however, about further classification of models with 
finitely many absorbing states. 

Models with infinitely many absorbing states constitute an important class. Such mod- 
els [§, ^, |ll], |I^ arises mainly in the study of surface catalysis but recently these models 
were related also with the Self-Organized Criticality |jl^ or biological evolution |]l4[. The 
critical behaviour of these models is very interesting. In the one-dimensional case, steady- 
state exponents have (1+1) DP values but certain dynamical exponents remain nonuniversal. 
However, some scaling arguments supported by numerical results suggest that the sum of 
these dynamical exponents is the same as in the (1+1) DP case [|l^]. 

In the two-dimensional case, fewer results are available and they are less accurate. For 
example, for the so-called dimer-dimer model Albano's estimations of the exponent (3 ~ 0.5 



describing the behaviour of the order parameter are marginally consistent with /3pp = 0.592 
of the (2+1) DP pB]. Similar calculations for the dimer-trimer model M, O] and for a certain 



variant of a sandpile model |13] also suggest the (2+1) DP universality class. Moreover, there 
are some renormalization-group arguments that models with infinitely many absorbing states 
should belong to DP universality class, at least with respect to the steady-state properties [^. 
The result which does not seem consistent with the (2+1) DP, namely /? ~ 0.2 — 0.22, was 



reported by Yaldram et al. for a certain model of CO-NO catalytic reaction |1C]. However, 
it was suggested by Jensen that their calculations were not accurate enough and additional 
calculations have shown that also in this case the model exhibits the (2+1) DP critical 



behaviour |12]. Indeed, Jensen calculated a number of exponents and all of them were 
consistent with the DP universality class. However, to estimate the exponent (3 and show 
that it is also consistent with the DP universality class, he used the so-called dynamic Monte 
Carlo method combined with some scaling relations, and such an approach is fundamentally 
different from the steady-state calculations by Yaldram et al. To refute Yaldram et al.'s result 
we would have be sure that indeed both approaches yield the same results. Although the 
dynamic Monte Carlo method is very frequently used and provides one of the most accurate 
estimations of critical parameters, its extent and the very reason of applicability, especially 
for models with infinitely-many absorbing states, is not, in our opinion, firmly established 
and further examples either supporting or contradicting the usage of this method would be 
very desirable. 

In the present Rapid Communication we study the two-dimensional version of a certain 



model which recently was introduced in the context of biological evolution |14]. This model 
has infinitely many absorbing states and its critical behaviour in one dimension [^] is in 
agreement with other models of this kind. Namely the critical exponent (3 ~ 0.273, describ- 



ing the density of an active phase, is very close to its (1+1) DP counterpart P-j^p ~ 0.2765 |18] 
and also the sum of dynamical exponents ij and 6 is universal (with respect to the choice 
of an absorbing state) and close to the DP value. We show, however, that in two dimen- 
sions the model has a number of unexpected features. First, our estimation of the exponent 
P ~ 0.273(5) is clearly different than its (2+1) DP counterpart /?^p ~ 0.592 |l|l. This value 



strongly suggests that in the steady state, due to a rather puzzling dimensional reduction, 
the critical behaviour of the two-dimensional model is the same as of its one-dimensional 



analogue |20|. Moreover, we show that the applicability of the dynamic Monte Carlo method 
to this model is highly questionable. In particular, we show that there exists a large class of 
absorbing states for which this method reproduces spurious dynamical transitions. We ex- 
pect that some other models with infinitely many absorbing states might also exhibit similar 
behaviour. Indeed, it was shown by Dickman that for certain two-dimensional model the esti- 
mation of critical point using the dynamical Monte Carlo method depends on the choice of the 



initial configuration and, in general, is different from the steady-state estimation [21]. How- 



ever, the reported difference was quite small (at most 4%) and analysis of critical behaviour 



is presumably affected by certain strong crossover effects j22|. Similar small differences were 
reported for yet another model with infinitely many absorbing states [^]. In the reported 
below results the difference between estimations of critical point using the dynamical Monte 
Carlo method and steady-state method might be as large as 30%. 

Our model is defined on a two-dimensional Cartesian lattice. Omitting the biological 
interpretation, we assign a certain number Wij to the bond connecting sites i and j and such 
that < Wij < 1. Introducing certain parameter r, we define the dynamics of our model as 
follows 1 14]: (i) Choose a site i at random, (ii) Calculate u = J2jWij, where summation is 



over all nearest neighbours j of the site i. (iii) If u; > r, then the chosen site i is active and 
all bond variables Wij are replaced by the new ones, chosen randomly. If a; < r, the chosen 
site is nonactive and its bond variables remain unchanged. It is obvious from the above rules 



that the model possesses infinitely many absorbing states. 

Since a computational implementation of the above rules is straightforward, we present 
only the results of our calculations. First we measured the density p of active sites (i.e., 
those with u; > r) in the steady state. The initial configuration of bonds is chosen randomly. 
Our results for various system sizes L are shown in Fig. Ij-Fig. |2|. These results clearly 
indicate the phase transition separating the active (p > 0) and absorbing (p = 0) phases of 
the model. Assuming that in the vicinity of the transition the density p has a power-law 
singularity p '-^ {vc — r)" and using the least-square method, we estimate re = 1.3867(5) and 
[5 = 0.273(5). These estimates are based on the results for L = 200 but the estimation of (5 
based on results for L = 100 is very similar. For L = 200 we made runs of 10^ Monte Carlo 
steps neglecting for each r data from the initial 10^ Monte Carlo steps. A Monte Carlo step 
is defined in a standard way, namely as a single, on average, update per site. 

Our results show that as far as the steady-state properties are concerned, this model does 
not belong to the (2+1) DP universality class. They instead strongly suggest that both one- 
and two-dimensional versions have the same exponents (3 as the (1+1) DP. At present we do 
not understand why such a dimensional reduction takes place. 

In our opinion, Yaldram et al.'s model might have the same critical exponent (3 as our 
model (i.e., we suggest that their calculations were inaccurate but not as much as suggested 



by Jensen). Why thus do Jensen's dynamic Monte Carlo calculations [12| yield the (2+1) DP 
behaviour? Although at present we cannot locate the cause of these inconsistencies, below we 
show that the dynamic Monte Carlo method requires serious reconsiderations when applied 
to models with infinitely many absorbing states. 

The idea of the dynamic Monte Carlo method is to set initially the model in one of the 
absorbing states with a seed of the active phase in the centre of the system and study the 
subsequent spreading of activity. One expects that for a certain value of the control parameter 
of the model various characteristics of spreading will exhibit a power-law scaling. Moreover, 
there is a considerable numerical support, mainly from studying one-dimensional models, that 
such dynamical critical point coincides with the steady-state critical point independently on 
the choice of an absorbing state. We show, however, that for two-dimensional models this is 
not the case and the choice of the absorbing state strongly affects the location and nature of 
the dynamical transition in the model. 

First, let us consider a trivial example, where as an absorbing state we have chosen a 
state with Wij = wq = Q for all bonds except the bonds surrounding a certain site, which 
are chosen such that this site is active. It is easy to realize that for r > 1 the activity 
cannot spread beyond this single site and the system quickly returns to the absorbing state. 
We do not present numerical data but we have checked that for r < 1 the activity usually 
spreads throughout the whole system, which indicates that the system is in the active phase. 
With such a choice of an absorbing state, r = 1 is the point which separates the active 
and absorbing regimes of the model. Although trivial, this is an example of an absorbing 
state for which the dynamical transition (r = 1) does not coincide with the steady-state one 
(r = 1.3867). It is also easy to notice that any absorbing state with i(;o € (0, ^^3-^) also yields 
spurious dynamical transition basically due to the same mechanism. 

It is also natural to expect that such spurious transitions will appear even when some 
inhomogeneous absorbing states are considered. As an example let us consider an absorbing 
state where all bonds Wi^j are chosen randomly from the interval (0,0.1). Using our previous 
analyses for homogeneous absorbing states one can see that for such a choice of an absorbing 



state the dynamic transition must take place at r^ € (1, 1.3) and numerical simulations |24] 

show that indeed r^ = 1.21(1), which is well below the steady-state critical point r^ = 1.3867. 

In the above examples the bond variables were set to low values and thus the active phase 

was strongly suppressed. Setting bond variables to large values (but such that their sums 



for each site do not exceed the threshold r), we can construct absorbing states where the 
absorbing phase is suppressed. As an example let us consider the case of wq = 0.25. Setting 
a central site in the active state, we measured the number of active sites N{t) as a function 
of time t and the results in the logarithmic scale for various r are shown in Fig. pi (average 
is taken over all runs). The number of independent runs varied from 200 for r = 1.42 to 
10000 for r = 1.47. It is essential in this type of simulations to ensure that the propagating 
activity never reaches the border of the lattice; for example for r = 1.45 we have to use 
L = 1500. It is also essential to keep a list of active sites since they constitute only a small 
fraction of all sites. One can see that r = 1.45(1) is a point which separates two regimes with 
asymptotically increasing and decreasing number of active sites. Why for Tc < r < 1.45, i.e., 
in the absorbing phase, might the activity spread for the infinitely long time? (Of course the 
activity can spread for the infinitely long time also for r < r^ but that is justifiable since for 
such r the system is in the active phase) . The reason for that is the unstable character of the 
absorbing state: large values of bond variables considerably ease spreading of activity. An 
example of such a propagating structure for r = 1.42 is shown in Fig. |^. We put a single site 
in the centre of the 1000x1000 lattice in the active state and recorded the configuration after 
time t = 1000. One can see that the activity is restricted only to the gradually increasing 
boundary of such a structure. But this is not surprising: for r = 1.42 the model is in the 
absorbing phase and activity in the centre dies out after some transient time which is needed 
for the system to find a stable absorbing state. One can also say that once such a structure 
has spread to infinity, an unstable absorbing state is transformed into the stable one. Let us 
also notice that the asymptotic slopes in Fig. y for r < 1.45 seem to be the same and slightly 
larger than unity. It suggests that structures like that shown in Fig. ^ might be fractals with 
the fractal dimension greater than unity and that this fractal dimension might be universal 
(i.e., independent of r). Similar propagating structures were observed also for other models 
with absorbing states pi] , 22]. 



We expect that wq = 0.25 is not the only value for which absorbing phase is suppressed 
and spurious dynamical transition is obtained. Similar results should be obtained also when 
a large-u; absorbing state contains some inhomogeneity (i.e., bonds are random variables 
from a certain range). But we also expect that there exists a large class of absorbing states, 
which are in some sense in between these two extremal classes examined above and for which 
dynamic Monte Carlo method will correctly locate the critical point. However, it means that 
in this method the choice of the absorbing state is very important and presumably it is very 
difficult to predict whether a given absorbing state will lead to a spurious or true critical 
point. 



In his calculations Jensen |12| used so-called typical absorbing states, which most likely 
correspond to our in-between states and which most likely correctly reproduce the transition 
point and yield the (2+1) DP exponents. We can only suggest that it is these extremal 
absorbing states which affect the steady-state dynamics and are responsible for the change 
of the universality class of our model. 
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Figure captions 

Fig. 1: Density of active sites p as a function of r calculated for L = 50 (A), 100 (o) and 200 

(+)• 

Fig. 2: The logarithmic plot of density p as a function of (re — r) for L = 200 and r^ = 1.3867. 

Fig. 3: The logarithmic plot of the number of active sites N{t) as a function of time t for (from 
the top to the bottom) r = 1.42, 1.43, 1.44, 1.45, 1.46 and 1.47. 

Fig. 4: Active sites (dots) propagating in the wq = 0.25 absorbing state for r = 1.42. A single 
active site was placed in the middle of the 1000x1000 lattice and the configuration 
shown was recorded after time t = 1000. 
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